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The molecules in lyotropic membranes are typically aligned with the surface normal. 
When these molecules are chiral, there is a tendency for the molecular direction to twist. 
These competing effects can reach a compromise by producing helicoidal defects in the 
membranes. Unlike thermotropic smectics, the centers of these defects are hollow and 
thus their energy cost comes from the line energy of an exposed lamellar surface. We de- 
scribe both the twist-grain-boundary phase of chiral lamellar phases as well as the isolated 
helicoidal defects. 
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1. Introduction 

Lyotropic liquid crystals share a number of features with their thermotropic cousins. 
In lyotropic lamellar systems, however, it is believed that because the molecules are am- 
phiphillic the molecular axes will align with the layer normal in an untilted, L a phase. Chi- 
ral molecules are well known to exert torques on each other [1] leading to chiral mesophases 
in thermotropic liquid crystals. In chiral lyotropics then, there will be a frustration between 
normal alignment and the tendency of the molecules to twist [2]. 

In thermotropic smectic-A phases, this frustration can be relieved via the formation 
of a twist-grain-boundary phase (TGB) [3] analogous to the Abrikosov flux line lattice of 
type-II superconductors. In this case the defects are screw dislocations with cores made of 
the associated nematic liquid crystal phase. In this letter we will describe a similar screw 
dislocation structure for lyotropic lamellae and propose an L a TGB phase. We will also 
discuss a possible defect mediated phase transition [4] between the L a phase and a normal 
cholesteric phase [5]. 

2. Helicoidal Defects in Membranes 

We model the free energy of an isolated bilayer membrane as a sum of contributions. 
The free energy for membrane fluctuations is [6] 



where R\ and R2 are the principal curvatures, N is the layer normal and V • N is the 
mean curvature. The free energy for fluctuations of the molecular axis n is given by the 
membrane version of the chiral Frank free energy, namely: 



where Vj_ refers to a gradient in the tangent plane of the membrane. At quadratic order, 
the non-linear terms arising from properly constructed, covariant derivatives do not come 
in. We couple these two free energies together by adding a term which favors the alignment 
of the layer normal with the molecular axis: 




(2.1) 




(2.2) 




(2.3) 
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The total free energy is the sum F = F m + F* + F mn . In the Monge gauge, we may 
represent the surface as a height function h{x, y) over the coordinates x and y. In this case 
the layer normal is N = [—d x h, —d y h, 1] / a/1 + (V±h) 2 . Expanding the nematic director 
about the ground state direction n = (z + 8n) / VT+ <5n 2 , the free energy to quadratic 
order in the deviations from the ground state becomes 

F = \j dxdy {7 (V±h + 5nf + k (Vi/i) 2 } + F*. (2.4) 

Note that if we require Sri = —V±h (which follows when 7 — > 00), or, in other words, if 
we require n||N, the resulting effective free energy has no chiral term since Vj_ x 5n = 0. 
It is only by relaxing the constraint that n be parallel to N that we see manifestations of 
chirality. Thus, free energies of the L a phase of membranes (with n||N) based only on the 
local membrane geometry will never show the effects of chiral molecules. The L a phase is 
the analog of the smectic- A phase in which all twist and bend have been expelled [7] . 

In the Lpi phases, there is an equilibrium projection c of the the director n onto 
the tangent plane, and there are manifestations of chirality in the free energy involving c 
[8,9,10], for example, through a term of the form J c • x c. These terms are responsible 
for a number of interesting modulated structures in and of the membrane. 

The free energy (2.4) is very similar to that for smectic-A liquid crystals, although 
the functions are restricted to depend only on the membrane coordinates. If we consider 
a screw-like configuration of the surface h(x,y) = 6arctan(y/a;) (where b is the "Burgers 
vector" of the defect), then we can calculate the energy of the associated configuration, 
after relaxing 8n0 = V ±h (i.e., n = N). This configuration, however, is singular at the 
origin and is not allowed. To facilitate the screw-like defect the membrane must cease 
to exist in the core, or, in other words, there will be exposed membrane edges around a 
solvent-filled cylindrical core. Thus, in the presence of a screw defect we must now add a 
line energy. This line energy comes from two pieces: the first is the energy cost of having 
an exposed edge, while the second comes from any surface terms in the free energy, for 
instance the Gaussian curvature. We group these energies together into a term: 

F e = n f d£, (2.5) 
Jc 

where C is the inner edge of the helicoidal surface. 

Far away from the cylindrical core, 5n = b [y, — x, 0] /(x 2 + y 2 ). Between the core 
radius £ and the twist penetration depth A = a/ -K2/7, the nematic director relaxes to its 
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large radius value. This means that the core has a double-twist texture, as found in liquid 
crystalline blue phases. A perfect helicoidal surface has no mean curvature, thus our defect 
does not contribute to F m . The free energies per turn are (with A = ^jK 2 /i): 



1 2 

mn = ~lb 7V In 

F n = -K 2 q b 



' A 2 ' 

1 + e. 



(2.6) 



F e = fXy/b 2 + (27T£) 2 , 



where, in the London (type II) limit, ^ <C A. To minimize the sum, note that only F n 
can be negative: this chooses the sign of b. Minimizing the sum of the constituent free 
energies, we find that there will always be a non-zero value of b for which the energy is 
In the two limits (6 -C £ or b 3> £) we find: 

27t£K 2 qo ,* f. 

27r 2 £ 7 ln(A 2 /£ 2 + l)+^ <<4 
b* = ( (2.7) 

U 7 ln(A 2 /£ 2 + l) >> ^ 
At the same time, £ is determined also via minimizing the free energy. We thus have 

7 2 A 4 (6*) 4 [(b*) 2 + (2nC) 2 ] = 16ttV [A 2 + (£*) 2 ] 2 (£*) 4 - ( 2 -8) 

The optimal value of £ = £* will determine the size of the hole in the membrane. 

Unlike a layered system where b is quantized in units of the spacing, b varies con- 
tinuously with K 2 qo- When a defect appears, the membrane must continue wrapping 
around into a many-layered helicoidal surface. Not to do so would create a line of exposed 
molecules leading out from the center of the defect r = to the boundary r = oo. 

The line energy \i presumably contains both energetic and entropic factors. In princi- 
ple, \i can be negative at the temperatures we are considering. If that were the case, even 
the non-chiral membrane, made of racemic mixture of molecules would unbind, leading to a 
randomly defective, layered structure. A lyotropic with a negative line tension would never 
really form a single layer structure and would have defects proliferating at the molecular 
level leading to a random unstructured collection of amphiphillic molecules. 

In the absence of any interactions that keep the membrane away from itself [11], there 
is no lowest energy state since the defect strength b need not be quantized for a single 
membrane. Entropic repulsion will balance against the energy gain of the defect, setting 
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the value of b. The entropic repulsion will scale as b~ 2 leading to some b* that will minimize 
the total free energy above the lower critical chirality. However, the entropic interaction 
will scale as the area of the bilayer whereas the line energy will scale like the length. The 
detailed balance between entropic repulsion and defect energy will thus be size dependent. 
It is also possible that a preferred, spontaneous torsion along the free edge would lead to a 
specific value for b [12]. Thus it should be possible for chiral lyotropic lamellae to exhibit 
isolated helicoidal fragments. 

3. Lamellar Melting and TGB Phases 

When membranes are stacked together, they can form multi-layered, lamellar struc- 
tures, similar to smectic-A phases in thermotropic liquid crystals. In highly swollen sys- 
tems, the layer spacing d is determined by entropic repulsion [11]. As in thermotropic 
smectics, the L a phase excludes twist [7] . Likewise, a sufficiently strong chirality can favor 
the entry of twist into the L a phase in the form of a TGB phase consisting of period- 
ically repeated twist-grain boundaries composed of periodic arrays of screw dislocations 
[3] . Though we know of no report of an experimental observation of this phase, we see no 
reason why it could not exist. 

Transitions in lyotropic systems are driven predominantly via changes in concentration 
rather than changes in temperature. We can determine the surfactant concentration, or 
equivalently the layer spacing d, at which the L a phase first becomes unstable with respect 
to a proliferation of dislocations by calculating the point at which the total energy per unit 
length of dislocation first becomes negative. This calculation is analogous to the calculation 
of Hci in a superconductor - it establishes the mean-field instability of the L a phase to 
the formation of a TGB phase. In the limit d ^> £, we find from (2.6) the total free energy 
per unit length of a dislocation is 



where the pre- factor of 1/d converts from free energy per turn to free energy per unit 
length [13]. 

Note that in (2.2), the elastic constants are membrane elastic constants, related to 
the bulk elastic constants by a factor of 1/d, the inverse layer spacing, and consequentially 
they do not change as the density changes (assuming that density changes are made most 




(3.1) 
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easily by changes in the layer spacing). The density is inversely proportional to the layer 
spacing, and so, against simple intuition, the energy of the defect (with b = d) decreases 
with increasing density. This is a simple result: as the layer spacing increases, the elastic 
energy of each distorted lamella grows since it must be distorted more to connect the 
consecutive layers. We thus predict that with increasing density the lamellar structure will 
be pocked with defects above a lower critical chirality. For a given value of K 2 qo — \i > 
there will be a critical spacing d* (critical density p* oc l/d*) below which (above which) 
the defects will penetrate. From (3.1), we find: 

7r 7 ln(A/0 

In general, a transition will occur when d is reduced to twice the preferred Burger's vector 
of the free membrane (2.7). If fluctuations are unimportant, then d* is the layer spacing 
at which there is a second-order transition from the L a to the TGB phase. 

Of course, it is possible that fluctuations are strong enough to destroy any TGB phase 
that might form and that the L Q phase transforms to a chiral N£ phase with twisted 
orientational order like that of a cholesteric. This phase is the liquid crystal analog of 
the melted vortex lattice of superconductors in a magnetic field [14] , which may intervene 
between the Meissner phase and the vortex phase of a type-II superconductor. It is a 
cholesteric phase formed by melting the TGB dislocation lattice while still retaining short- 
range smectic order. Heat capacity measurements provide strong evidence of this phase 
in thermotropic systems [15]. Possible phase diagrams based on the analogy between this 
system and high-Tc superconductors are shown in Figure 2. 

We propose that the layers melt via dislocation loop unbinding, as in the smectic- 
A-to-nematic transition [4]. Here, however, the chiral bias will cause one handedness 
of screw dislocation to be preferred over the other. This bias should change a second- 
order-like unbinding transition to a first-order transition: microscopic defect loops can 
no longer unbind smoothly since a loop must contain an equal number of left and right 
handed screws. Either the defect loops will unbind and the dissident part of the loop will 
move to the boundary, or, equivalently, the appropriate dislocations will nucleate from the 
boundary. 

The chiral interaction also changes the simple unbinding picture in another way: from 
the above we can see that the free energy difference per unit length between a wrongly- 
handed screw and a correctly- handed screw will be: 

AF/L = 2\K 2 k b\. (3.3) 
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We thus expect that, in the presence of thermal fluctuations, at the unbinding transition 
of the defect loops, the largest size dislocation loop will have a characteristic length 

L " = w- (3 ' 4) 

As the chirality of the molecules becomes small, i.e. as k — > 0, the screw dislocations will 
grow to infinite size and will proliferate, as in the usual non-chiral scenario. 

We can estimate the latent heat per molecule via a scaling argument. If we assume 
that without chirality there is a second-order smectic-A-to-nematic transition (or nearly 
second order), then the transition will be rounded out when the correlation length is on the 
order of L* . In this fluctuation dominated "type II" limit the transition from smectic order 
to isotropic order would proceed via an intermediate cholesteric-like phase (either TGB or 
iV£) as shown in Figure 2(a). We believe that the most likely candidate for this phase is 
the iV£ phase [15] of thermotropic smectics [16]. It is also possible that the isotropic phase 
intervenes before a TGB phase ever forms, as is shown in the phase diagram of Figure 
2(b). For fixed d* we may calculate A 2 /£ 2 . In the limit that d* > and A > £*, we have: 

^ 1 = ^L_ 3 . (3.5) 

Thus we expect extreme type II behavior, i.e. large when -ft^Qo ~ A* or when -KT27 is very 
large. The transition from the lamellar to TGB phases would presumably be second-order 
while the transition from the lamellar to the iV£ phase would be first-order due to the 
strong fluctuations. 

In the latter case we can estimate the scaling behavior of the free energy per unit 
volume at the transition by cutting off a the transition at the size L*. Scaling gives us for 
the entropy per unit volume: 

df r 



(L*)" 1 /" ~ \K 2 k \ {1 - a)/ \ (3.6) 



where / is the free energy density, t is the reduced temperature, v is the correlation length 
exponent and a is the usual specific heat exponent of the nearby second-order transition. 
Note that this result based on dislocation loop unbinding gives the same scaling result in 
terms of ko as that obtained previously via Landau theory [17]. 

Since the defects in the lamellar phase are necessarily accompanied by holes in the lay- 
ers, X-ray scattering should provide a view of the dislocation unbinding transition leading 
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to a chiral TGB or cholesteric phase. As the dislocation loops grow to L* , the scattering 
will show liquid like peaks at q± = 47r/(y / 3L*) oc kod, where q± is the wavevector perpen- 
dicular to the layer normals. Though there will be a proliferation of holes in each layer, 
only the holes coming from the same dislocation loop will be correlated, leading to an 
X-ray structure with peaks at the characteristic inverse length scale. This could be tested 
by changing the concentration of chiral molecules (or, equivalently, adjusting the ratio of 
left-handed to right-handed enantiomers), as well as by adjusting the layer spacing. 

Finally we mention that it is possible for the lamellar stacks to be in a "type I" limit. 
In this case the melting from the stacked phase to the isotropic phase will occur at one 
place: there will be no intervening nematic or cholesteric-like phase. We would expect this 
transition to be first-order and would not occur via dislocation loop unbinding. Instead we 
would expect this transition to occur via dislocations nucleating at the boundary. We also 
note that in the zero chirality limit there are no known transitions from lamellar structures 
to nematic structures. In the chiral case, however, we have seen that chiral mesophases 
arise via the presence of chiral topological defects. 
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Figure Captions 



Fig. 1. A screw dislocation in a lyotropic lamellar phase. Note that the core is devoid 
of all molecules and is filled with solvent. In thermotropic smectics, the core is 
nematic and contains the mesogens. 

Fig. 2. Possible phase diagrams for the lyotropic system as a function of chirality qo 
and surfactant density p. Solid lines indicate first-order transitions while dashed 
lines indicate second-order (or weakly first-order) transitions, (a) Lyotropic with 
four distinct phases: lamellar, TGB, iV£ and isotropic, (b) It is possible that 
the isotropic phase intervenes before the TGB phase occurs. The phantom TGB 
region is shown with hatched lines: it will never appear. 
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